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Abstract. Fully terrestrial lifestyles have evolved in a lineage of salamanders living since the Oligocene in the Western 
Palaearctic that adaptively radiated in morphology, ecology and reproductive modes. I analyse the relationships between 
body size and shape, and climate and reproductive mode in true terrestrial salamanders of the family Salamandridae. 
Principal component analyses were carried out on bioclimatic and morphologic data and used to reconstruct the diver- 
sification of these along the evolutionary timeline of this lineage. Additionally, a phylogenetic generalized least squares 
regression was performed in order to test the effect of climate and reproductive mode on their morphologic diversifica- 
tion. Western Palearctic salamanders experienced a strong disparity between small and elongated forms early during their 
evolution, and short-bodied species exhibit large size variation. Body elongation in the Salamandridae seems to be linked 
to an adaptation to semifossorial habits and is currently restricted to a narrow range of climatic conditions, which is in 
contrast to the other terrestrial species. Within this latter group, large species exclusively live at the southernmost edge of 
the geographic range of true terrestrial salamanders. Thus, Bergmann’s rule is not adhered to by true salamanders, since 
larger species live in hot climates. The ancestral oviparous mode of reproduction in the Salamandridae could be main- 
tained in elongated salamanders due to their body shapes not posing constraints, whereas the largest sized short-bodied 


species have developed larviparism. 
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Introduction 


Based on molecular evidence and fossil records, Salaman- 
drids appeared 113-55 Mya (ESTES 1981, ZHANG et al. 2008) 
and diversified into 105 species distributed in North Amer- 
ica, Europe and temperate, but also tropical, regions of 
Asia. The ecological radiation of this family involved di- 
versifications of body shape, tongue morphology, repro- 
ductive strategies, and courtship behaviour (e.g., OZETI & 
WAKE 1969, WAKE & OZETI1969, HOUCK & ARNOLD 2003). 
Remarkably, one salamandrid lineage has developed sev- 
eral morphological and reproductive innovations in order 
to shed its dependence on water, leading to truly terrestrial 
lifestyles in some species (MILLER & LARSEN 1990, WAKE 
2004). Named “true” salamanders (VEITH et al. 1998), they 
live in the western parts of the Palaearctic and are repre- 
sented by the genera Chioglossa, Mertensiella, Salamandra 
and. Lyciasalamandra. Phylogenetic evidence suggests an 
Oligocene origin of the Western Palaearctic salamanders 
that could be placed in the Near East (RopDRIGUEZz et al. 
2017), where Mertensiella, Lyciasalamandra, and one basal 
lineage of Salamandra, S. infraimmaculata live (Ent et al. 
2019). Since that period, they have experienced diversifica- 
tions of body size, shape and coloration, colonizing a wide 


variety of climates ranging from that of Mediterranean is- 
lands to that of Alpine mountain ranges. However, only the 
species of Lyciasalamandra have become totally decoupled 
from water by developing not only terrestrial reproduc- 
tion, but also viviparity, even supporting it with matro- 
trophy (HAFELI 1971, BUCKLEY et al. 2009). Accordingly, it 
is likely that in Western Palaearctic salamanders, ecology, 
body size and shape evolved driven by the disparity of cli- 
matic environments and the reproductive strategies devel- 
oping in this lineage. 

Correlations between the evolution of body shape and 
other biological traits, or even abiotic ones, has been re- 
ported by several researchers and is considered evidence 
of adaptation (e.g., GRIZANTE et al. 2012, BONETT & BLAIR 
2017). The main aim of this study is to characterize the ex- 
tent of body size and shape diversity in true salamanders, 
how it evolved through time, and test whether current cli- 
matic conditions and reproductive modes may have acted 
as drivers of their morphologic diversification. 


Material and methods 


A sample of 768 specimens of Western Palaearctic salaman- 
ders, comprising all species of Mertensiella, Salamandra, 
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Chioglossa, and most of Lyciasalamandra (Fig. 1) from col- 
lections housed at natural history museums was analyzed 
(Supplementary document S1). In order to capture the vari- 
ation of body size and shape, seven variables were taken: 
length from the tip of the snout to the posterior margin of 
the cloacal protuberance (SVL), head length and maximum 
width, forelimb, hindlimb and tail length, and inter-limb 
distance measured between the closest point of insertion of 
extremities on the right ventral side (Table 1). Despite the 
sexual dimorphism in some traits (REINHARD et al. 2015), 
differences between lineages of terrestrial salamanders are 
likely more phylogenetically informative than those be- 
tween sexes. Therefore and in order to maximize sample 
sizes, males and females were pooled for the analysis. Bio- 
climatic data on each taxon were extracted at 30 arc-second 
intervals (7182 records) from the Worldclimate database 
(Fick & HIjMANS 2017), using gathered data on their dis- 
tribution (Herpnet, Inaturalist, and other sources: Supple- 
mentary document S2). The following variables represent- 
ing climatic trends and extreme variations were selected: 
mean annual temperature, maximum temperature of warm- 
est month, minimum temperature of coldest month, annual 
precipitation, precipitation of wettest quarter, and precipita- 
tion of driest quarter. Principal component analyses (PCAs) 
of morphometric and bioclimatic data sets were performed 
using the mean scores for each salamander taxon with the 
aid of Statistica 4.5. In order to test whether morphological 
and climatic variables are labile or strongly shaped by the 
evolutionary history, the phylogenetic signal was checked 


using the K and Lambda tests (PAGEL 1999, BLOMBERG et al. 
2003) incorporated in the function phylosig in Phytools 0.6 
(REVELL 2012) using R 3.6.2. A phylogenetic tree of Western 
Palaearctic terrestrial salamanders was built using available 
information about their evolutionary relationships (Ro- 
DRIGUEZ et al. 2017, EHL et al. 2019). Branch lengths were 
estimated by mean maximum likelihoods upon the HKY- 
86 model with invariant sites and gamma modelling of 
variable sites using a dataset formed by 12sRNA and cyto- 
chrome b mitochondrial genes (Supplementary document 
S3). Branch lengths were ultrametrized by applying a mo- 
lecular clock based on three external points of calibration 
(Supplementary document S4) applying the RelTime meth- 
od (TaMuRA et al. 2012). All these analyses were performed 
by means of Mega 7.0 (KUMAR et al. 2016). 

Evolution of morphologic and climatic PCA was traced 
by means of character reconstructions based on maximum 
likelihoods using the function fastAnc incorporated in 
Phytools 0.6 and the ultrametric phylogenetic tree. To test 
whether climate and reproductive mode affected body size 
and shape diversification, a phylogenetic generalized least 
squares (PGLS) regression was performed. Reproductive 
mode was defined by means of three categorical variables, 
oviparity, larviparity and pueriparity, coding the absence 
or presence of a trait as o or 1, which facilitated the inclu- 
sion of taxa with interpopulation variation (i.e., Salaman- 
dra salamandra gallaica, S. s. bernardezi). Analyses were 
conducted using Phytools 0.6 and non-phylogenetic line- 
ar-squared regressions on R 3.6.2. 


Figure 1. Geographic distribution of the four genera of Western Palaearctic salamanders: monotypic and elongated Chioglossa (salmon) 
and Mertensiella (orange), and short-bodied Lyciasalamandra (black) and Salamandra (blue). The latter genus is geographically the 
most widespread of the true salamanders. 
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Table 1. Means and standard deviations of morphometric variables of Mertensiella, Chioglossa, Lyciasalamandra and Salamandra. 


Taxa N SVL Head Head Limb Forelimb Hindlimb Tail 
length width interval length length length 
Chioglosa lusitanica 5 48.5+0.7 9.4+0.3 5.6+0.4 25.041.6 9.545.2 12.6+0.4 86.0+8.1 
Mertensiella caucasica 62.545.3 14.4+1.1 9.1+0.6 30.744.0 16.3+0.7 19.641.0  84.7+14.2 
Lyciasalamandra fazilae 10  60.07.9 14.0+1.4 10.4+1.1 27.0+4.2 17.8+2.7 20.6+2.3 47.7+7.7 
Lyciasalamandra luschani 10 66.944.0 14.9+0.9 10.8+0.7 31.4+4.4 19.1+1.1 22.041.6  48.644.3 
Lyciasalamandra atifi 2 74.9+13.3 17.3+1.6 12.041.9 35.446.9 21.642.7 25.242.1  64.5+12.0 
Lyciasalamandra billae 2- 67.941.5 14.9+0.1 10.6+0.2 30.741.9 17.641.2 21.7+2.0 53.442.7 
Lyciasalamandra flavimembris 2  66.740.4 15.0+0.4 10.6+0.2 31.7+40.6 18.8+3.1 21.5+0.7 52.9+3.3 
Salamandra atra 27 76.3+15.2 17.343.1 12.042.4 35.24+8.9 20.4+4.9 23.6+4.5 50.4+11.0 
Salamandra corsica 18 100.8+8.3 23.3421 18.3+1.4 44.0+4.9 29.4+3.8 32:9%3.9 60.0+8.2 
Salamandra i. infraimmaculata 8 128.2+22.2 26.644.0 21.3+3.1 60.5+12.0 33.4+6.2 38.147.9  82.8+14.8 
Salamandra lanzai 9 = 83.7+5.0 20.4+1.0 14.3+0.7 37.642.9 21.641.3 25.5+1.4 58.3+6.1 
Salamandra a. algira 97.941.5 22.4+0.4 15.8+0.3 47.343.8 27.141.3 29.8+0.2 64.745.8 
Salamandra a. splendens 3  74.6+4.9 17.4+0.9 11.9+0.6 35.643.1 19.3+0.4 22.2+2.2 55.9+2.9 
Salamandra a. spelaea 15- 89.3416.5 = 20.342.6 14.5+1.8  42.1+10.1 23.7+4.7 27.5+45.5  63.2+17.4 
Salamandra a. tingitana 12 95.3+9.3 21.041.6 15.241.1 44.6+4.7 25.2+2.5 30.643.9 71.5+8.6 
Salamandra a. atra 27 76.3+15.8 17.343.1 12.0+2.5 35.049.1 20.4+5.0 23.644.6  50.711.1 
Salamandra s. almanzoris 75 77.0+7.1 18.4+2.1 14.8+1.5 33.2+4.5 21.4+2.7 24.4+2.7 49.4+6.0 
Salamandra s. bejarae 48 97.1+8.2 21.8+2.2 16.741.5 43.345.4 22.743.1 30.942.9 59.7+8.2 
Salamandra s. bernardezi 27  75.4+9.7 16.6+1.7 12.1+1.6 33.145.0 19.2+2.8 22.9+2.7 47.8+6.6 
Salamandra s. crespoi 2  91.6+6.0 22.0+0.9 15.7+£0.7 42.0+5.0 24.8+3.1 30.9+0.9 55.643.6 
Salamandra s. fastuosa 113 82.6+10.3 18.1+2.4 13.741.8 35.3+6.1 22.4+3.2 26.0+3.7 52.4+8.3 
Salamandra s. gallaica 81 98.5+10.1 22.4+2.1 17.3+1.8 42.7+6.1 28.0+3.6 32.143.8 59.7+8.1 
Salamandra s. giglioli 40 85.9+6.4 19.6+1.6 14.841.2 373E39 24.1+2.4 27.242.6 54.6+6.4 
Salamandra s. longirostris 36 110.2+9.7 24.2+2.6 18.0+1.5 51.9467  -28.9+3.5 34.1+3.5 64.9+7.9 
Salamandra s. morrenica 52 103.0+48.9 23.2+2.0 17:0+1.9 47.2+4.9 28.8+3.6 32.9£3.9 61.7+8.7 
Salamandra s. salamandra 9 95.2107 2192.5 16.1+2.1 43.9+7.0 26.0+3.0 29.9+3.3 62.5+8.3 
Salamandra s. terrestris 120 95.1+9.6 21.2+2.0 15.8+1.6 43.2+6.2 25.642.6 30.043.2 58.7£7.5 
Results Table 2. PCA loadings of mean values in the seven linear-mor- 


Body size diversified during the evolution of terrestrial sal- 
amandrids and currently ranges from the smallest species 
Chioglossa lusitanica (maximum SVL 49.2 mm) to the larg- 
est, Salamandra infraimmaculata (max. SVL 150.4 mm, Ta- 
ble 1). Body size and shape diversifications were well rep- 
resented by the PCAs, accounting for a larger portion of 
variation of the seven morphometric variables (Table 2). 
The first PC accounted for 82.8% of the total variation and 
displayed a typical size structure in the shape of high and 
negative coefficients in all the morphometric variables. The 
second PC was polarized by the negative coefficients of tail 
length and inter-limb distance, and on the other hand by 
the extreme positive coefficients of head width and fore- 
limb length. Thus, this second component may be likely 
associated with body shape diversity. Western Palaearctic 
terrestrial salamanders displayed two contrasting patterns 
of body shape (Fig. 2a). Chioglossa and Mertensiella are two 
closely related small-bodied lineages characterized by slim 
and elongated bodies, tails and narrow heads, a trend more 
developed in the North Western Iberian than in the Cau- 
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phometric variables covering the morphological diversification 
of Western Palaearctic salamanders. 


Variable 1* PC 24 PC 
SVL -0.993 -0.059 
Head length -0.988 0.017 
Head width -0.982 0.137 
Limb interval -0.956 -0.163 
Forelimb length -0.978 0.136 
Hindlimb length -0.991 0.066 
Tail length -0.138 -0.987 
Eigenvalues 5.802 1.048 
% of explained variation 82.886 97.863 


casian endemics. Body shape in Lyciasalamandra and Sala- 
mandra was roughly similar, and the main trait explaining 
interspecific variation was overall SVL. Thus, Lyciasala- 
mandra was turned up as a genus of uniformly small-sized 
salamanders grouped together with Salamandra atra and 
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the subspecies bernardezi, giglioli and almanzoris of S. sala- 
mandra, all of which are small-sized as well. At the opposite 
end, the Middle Eastern Salamandra infraimmaculata and 
the North African S. algira have the largest SVLs, followed 
by the southern Iberian subspecies S. s. longirostris. Based 
on their placement in the second PC, Salamandra infraim- 
maculata exhibited a slight trend towards body elongation. 

Bioclimatic variables selected for the analyses reflected 
the diversity of environmental conditions experienced by 
Western Palaearctic salamanders. The first PC explained a 
total variation of 60.5% (Table 3), displaying an inverted in- 
fluence of temperature vs. precipitation on the first PC. The 
second PC produced high negative coefficients for annu- 
al precipitation and maximum precipitation in contrast to 
the positively weighted minimum precipitation. The harsh- 
est environmental conditions (i.e., lower precipitation and 
high temperatures) were found in the Middle East, which 
is inhabited by Lyciasalamandra and Salamandra infraim- 
maculata and, to a lesser extent, by Salamandra algira 
(Fig. 2b). Interestingly, the two elongated salamander gen- 
era Chioglossa and Mertensiella exhibited close affinities to 
climates characterized by buffered temperatures and high 
amounts of precipitation. At the other extreme end of the 
climatic variation experienced by true salamanders, cold 
and rainy climates were characteristic of the habitats of the 
Alpine species Salamandra lanzai and S. atra. 

Body size and shape diversifications in Mertensiella, 
Chioglossa, Lyciasalamandra and Salamandra were in- 
fluenced by the evolutionary history of this grouping, al- 
though the strength of the phylogenetic effect strongly dif- 
fers between methods of testing. The Lambda test suggest- 
ed strong phylogenetic constraints on the evolution of size 
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Table 3. PCA performed on the six selected bioclimatic variables 
along the distribution of Mertensiella, Chioglossa, Lyciasalaman- 


dra and Salamandra. 


Variable 1* PC 24 PC 
Mean annual temperature 0.964 -0.132 
Max temperature of warmest month 0.967 -0.015 
Min temperature of coldest month 0.875 -0.019 
Annual precipitation -0.459 -0.871 
Precipitation of wettest quarter 0.203 -0.964 
Precipitation of driest quarter -0.865 0.051 
Eigenvalues 3.633 1.710 
% of explained variation 60.566 89.076 


and shape (first PC: A=0.867, P<o.001; second PC: \=0.922, 
P<o.001) while the K test turned up a weak effect (first PC: 
K=0.014, P=0.028; second PC: K=o.029, P=o0.011). Similar 
results were obtained in the case of the first PC extracted 
from bioclimatic variables (A=0.966, P<o.001; K =0.015, 
P=0.037), but not in the second PC (A = 0.241, P=0.423; 
K=0.002, P=0.400). Therefore, only non-phylogenetic lin- 
ear-squared regressions were performed on this PC. 

The ancestral condition of terrestrial Western Palaearc- 
tic salamandrids was not that of an elongated body (Fig. 3a), 
which was first developed by the common ancestor of Chi- 
oglossa and Mertensiella. Based on maximum likelihood 
reconstruction, the common ancestor of this group was a 
small-sized and short-bodied species (Fig. 3b). From this 
ancestral condition, SVL slightly increased in the common 
ancestor of Lyciasalamandra and Salamandra. The former 
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Figure 2. Morphological and climatic diversification in terrestrial Western Palaearctic salamanders: A) Morphospace generated by 
the two first PCs extracted from seven linear-morphometric measurements. B) Range of variation of climatic conditions within the 
range of terrestrial Western Palaearctic salamanders as represented by the two first PCs extracted from the six bioclimatic variables. 
Taxa are symbolized as: full rhombs, Chioglossa lusitanica (blue), Mertensiella caucasica (red); asterisks, L. atifi (black), L. billae (red), 
L. fazilae (green), L. flavimembris (cyan), L. luschani (blue); empty triangles, Salamandra algira algira (red), S. a. splendens (black), 
S. a. tingitana (cyan), S. a . spelaea (green), Salamandra infraimmaculata (blue); full squares, Salamandra atra (red), S. lanzai (blue), 
S. corsica (black); full circles, Salamandra s. almanzoris (cyan), S. s. salamandra (red), S. s. terrestris (black), S. s. bejarae (green), S. s. 
gallaica (blue); empty squares, S. s. bernardezi (blue), S. s. fastuosa (red), S. s. giglioli (black); empty circles, S. s. crespoi (blue), S. s. 


longirostris (red), S. s. morenica (black). 
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Figure 3. Maximum likelihood reconstruction of the evolution of morphology based on the first PC (a) and second PC (b) in ter- 
restrial Western Palaearctic salamanders. 
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genus then retained this medium size while Salamandra 
experienced further body length diversification. Large 
salamanders independently appeared several times dur- 
ing the evolutionary history of this genus, i.e., in Salaman- 
dra corsica, S. s. longirostris, the crespoi-morenica lineage, 
and today’s largest species, S. infraimmaculata. Converse- 
ly, size reduction evolved in the Ibero-south Italian lineage 
of Salamandra salamandra, as represented today by S. s. 
almanzoris, giglioli and bernardezi, and the Alpine clade of 
S. atra and S. lanzai. Maximum likelihood reconstruction 
indicates that the original climatic conditions in which the 
common ancestor of Salamandra, Lyciasalamandra, Chio- 
glossa and Mertensiella evolved were similar to those cur- 
rently experienced by these two small and elongated lat- 
ter species (Fig. 4). Adaptation to xeric environments oc- 
curred independently five times during the evolution of 
the group (in Lyciasalamandra and four lineages of Sala- 
mandra) while extreme cold habitats were colonized only 
once during the evolution of the Alpine lineage. 

The evolution of body size and shape in Western Palae- 
arctic terrestrial salamanders was found to be related with 
the climatic conditions. PLGS regression revealed that first 
PC extracted from the six bioclimatic variables was sig- 
nificantly related with size and shape PCs (Table 4). The 
relationship was positive with the first morphologic PC, 
revealing that SVL increased with rising temperature and 
decreased with rising precipitation. The opposite pattern 


was found in the second PC, indicating that body elonga- 
tion was related with buffered precipitation and tempera- 
ture during the driest and coldest months. Linear-squared 
regression of morphologic PCs on the second climatic PC 
was not significant in any of the tested scenarios (Table 4), 
congruent with the low percentage of variation accounted 
for by this component. PGLS regression exhibited some re- 
productive modes to have a significant effect of on body 
size and shape (Table 5). Thus, larviparous salamanders 
tend to be largest, but the other modes did not show a sig- 
nificant relationship with body size, and oviparity is related 
with body elongation. 


Discussion 


The wide array of climatic conditions experienced by West- 
ern Palaearctic salamanders during their long evolutionary 
history (EHt et al. 2019) may have driven their diversifica- 
tion in body shape and size, as my results suggest. Early 
during their evolution, true salamanders strongly diverged 
into elongated and short-bodied forms probably as a re- 
sult of their ecological divergence in a quest to exploit dif- 
ferent niches. Chioglossa and Mertensiella thus developed 
semi-fossorial habits and now live in microhabitats with 
high moisture levels generated by the proximity of fast- 
flowing streams in mountain areas (TARKHNISHVILI 1993, 
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Figure 4. Maximum likelihood reconstruction of the evolution of first principal bioclimatic component among terrestrial Western 


Palaearctic salamanders. 
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Table 4. Relationships between bioclimatic (CPC) and morphometric (MPC) principal components of Western Palaearctic salaman- 
ders assessed by PLGS and OLS regressions. Significant p-values are given in bold. 


Predictor Predicted Slope + SE T P Intercept + SE T P R? 

is CpG 1* MPC 0.048+0.020 2.357 0.029 0.001+0.0020 0.001 0.999 0.235 
2™4 MPC -0.050+0.205 -2.478 0.023 0.001+0.019 0.001 1.000 0.254 

2” CPC 1* MPC -0.02540.019 -1.324 0.197 0.001+0.0189 0.001 1.000 0.065 
274 MPC 0.019+0.019 0.995 0.329  -0.001+0.019 0.001 1.000 0.038 


Table 5. Relationships between morphometric (MPC) principal components and reproductive mode of Western Palaearctic salaman- 


ders assessed by PLGS regression. 


Predicted Predictor Slope + SE T P Intercept + SE T P R? 
Oviparity 2.011+1.062 1.892 0.070 0.040+1.105 0.036 0.970 0.129 

1* MPC Larviparity -1.346+0.439 -3.061 0.005 1.321+0.839 1.574 0.128 0.280 
Puelparity 0.464+0.270 1.721 0.098 0.920+1.151 0.799 0.431 0.109 
Oviparity -2.610+1.007 -2.589 0.016 0.168+1.065 0.157 0.875 0.218 

2°? MPC Larviparity -0.583+0.559 -1.042 0.307 -1.087+1.283 -0.847 0.405 0.043 
Puelparity 0.4154+0.211 1.959 0.061 -1.348+1.308 -1.030 0.313 0.137 


ARNTZEN 1995). Body elongation could be an adaptation 
to this semi-fossorial lifestyle, as it is also found in other 
salamanders with similar lifestyles (for example, Batracho- 
seps or Oedipina, PARRA-OLEA & WAKE 2001), but at the 
same time it may impose strong physiological constraints. 
Elongated bodies have a greater body surface/volume ra- 
tio (BAKEN et al. 2019) that will increase the loss of water, 
which is dramatically important in species that breathe via 
their skin. Hence, elongated Salamandridae are inevita- 
bly restricted to very humid microhabitats. However, our 
analyses showed that the first climatic PC was influenced 
by the phylogeny and a strong climatic similarity between 
Chioglossa and Mertensiella. Fossil and biogeographic evi- 
dence indicates that the ancestor of these elongated sala- 
manders had a large geographic distribution in the past, 
becoming restricted to climatic conditions currently found 
only in northwestern Iberia and the south of Caucasus 
(BLAN et al. 2009, EHL et al. 2019). Hence, this climatic 
affinity between the two species could be the result of phy- 
logenetic conservatism (WIENS & GRAHAM 2005). In con- 
trast, Lyciasalamandra and Salamandra have shorter bod- 
ies and live on the floor, mainly in forested habitats. In or- 
der to avoid the risk of desiccation they have developed 
almost strictly nocturnal activity patterns during rainy or 
otherwise humid nights. 

Remarkably, my analyses revealed that Western Palae- 
arctic salamanders defy Bergmann’s rule (WATT et al. 2010) 
and exhibit an opposing trend. Very limited evidence for 
the validity of this rule has been found in amphibians and 
in salamanders in particular (ASHTON 2002, ADAMS et al. 
2008). Likely, true salamanders deviate from this rule be- 
cause the largest salamander species (Salamandra algira 
and S. infraimmaculata) happen to inhabit the warmer and 
xeric areas at the southernmost margins of the overall dis- 
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tribution range. Even though high-altitude taxa such as, for 
example, Salamandra atra and S. s. almanzoris, are smaller, 
and the largest member of the Alpine clade, S. lanzai, lives 
at the highest altitudes in the coldest climates (ANDREONE 
& SINDACO 1989). 

The harsh climatic conditions faced by salamanders at 
the southernmost edge of the lineage distribution could 
force them to strongly reduce their annual activity (DE- 
GANI 1996), thus producing delayed sexual maturity, con- 
sequently facilitating a longer period of subadult growth, 
and therefore leading to larger sizes. Nevertheless, the 
delay in reaching sexual maturity is not so large as to ex- 
plain the differences in SVL between large and small spe- 
cies (SINSCH et al. 2017, ALTUNISIK 2018), thus suggesting a 
more complex background. Perhaps large salamanders can 
adapt better to dry habitats by having a lower ratio between 
body surface and volume. 

Viviparity is one remarkable evolutionary process to 
adapt to the truly terrestrial lifestyle that is today found 
among urodeles only in Western Palaearctic salamanders 
(BUCKLEY et al. 2007). Within this group, oviparity only 
occurs in the elongated Mertensiella and Chioglossa spp., 
while the short-bodied Lyciasalamandra and Salamandra 
spp. have larviparous and pueriparous modes of repro- 
duction (BUCKLEY et al. 2009). Thus, it is likely that the 
development of small and elongated bodies in Chioglossa 
and Mertensiella has enforced the retention of the ancestral 
oviparous mode of reproduction by not providing enough 
abdominal space for bearing larvae. Despite this effect of 
body size on the trade-off between clutch and offspring size 
(SALTHE 1969), my results provide very limited support for 
a clear relationship with larviparous or pueriparous repro- 
duction. Species and subspecies at the southern edge of the 
geographic area inhabited by Salamandra are those with 
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the largest sizes and are mainly larviparous, but the north- 
ernmostly distributed Salamandra algira tingitana also 
gives birth to numerous fullly developed young (DONAIRE 
et al. 2001). Pueriparity is obligatory in the Mediterrane- 
an and small-sized Lyciasalamandra (OZETI 1979), but it is 
also the mode of reproduction of the Alpine salamanders 
and some populations of Salamandra Salamandra, namely 
the subspecies gallaica, bernardezi and fastuosa (GARCÍA- 
Paris et al. 2003, VELO-ANTON et al. 2012). The develop- 
ment of the variety of derived reproductive modes, from 
larviparity to pueriparity with oophagy, adelphophagy or 
matrotrophy, seem to constitute very complex phenomena 
driven by factors other than body size (Dopazo & KOREN- 
BLUM, 2000, VELO-ANTON et al. 2007, 2012, BEUKEMA, 
2012, ALARCON-Rios et al. 2020). 

As my results revealed, within the true salamanders, 
Salamandra salamandra has the most intraspecific variable 
body size, making it an excellent model species. The study 
of interpopulation differences of demography, reproduc- 
tive parameters, and body size in the light of their complex 
biogeography, might be helpful to understanding interrela- 
tionships between causative factors that shape body size in 
true salamanders. 
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